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ABSTRACT 

In this paper, we focus on statistical region-based active con- 
tour models where image features (e.g. intensity) are ran- 
dom variables whose distribution belongs to some paramet- 
ric family (e.g. exponential) rather than confining ourselves 
to the special Gaussian case. Using shape derivation tools, 
our effort focuses on constructing a general expression for 
the derivative of the energy (with respect to a domain) and 
derive the corresponding evolution speed. A general result is 
stated within the framework of multi-parameter exponential 
family. More particularly, when using Maximum Likelihood 
estimators, the evolution speed has a closed-form expression 
that depends simply on the probability density function, while 
complicating additive terms appear when using other estima- 
tors, e.g. moments method. Experimental results on both syn- 
thesized and real images demonstrate the applicability of our 
approach. 

1. INTRODUCTION 

In image segmentation, the main issue is to extract one or 
several regions according to a given criterion. Region based 
active contours have proven their efficiency for such a task. 
The evolution equation of the deformable curve is deduced 
from a functional to minimize that may take benefit of statis- 
tical properties of the image features. In this case, the func- 
tional depends on the probability density function (pdf) of the 
image feature within the region. For instance, the minimiza- 
tion of the - log-likelihood of the pdf has been widely used 
for the segmentation of homogeneous regions in noisy images 
lUElO. In these works, the pdf is parametric, i.e. it follows 
a prespecified law (Gaussian, Rayleigh). The pdf is then in- 
dexed by one or more parameters (e.g. mean and variance for 
Gaussian laws) describing the distribution model. These pa- 
rameters depend on the region and must be estimated at each 
evolution step. The estimation may be performed using differ- 
ent methods (moments method, maximum likelihood (ML)). 
However, to the best of our knowledge, the influence of the 
estimation method on the computation of the evolution equa- 
tion has never been investigated. 



In this paper, using shape derivative tools |l4l|5]|6l, our 
effort focuses on constructing a general expression for the 
derivative of a functional depending on parametric pdfs, and 
hence on the corresponding evolution equation. Our contri- 
bution is twofold. Firstly, we show that the estimator of the 
distribution parameters is crucial for the derived speed expres- 
sion. Secondly, we propose to give general results for the 
evolution equation within the framework of multi-parameter 
exponential family and ML estimation theory. Therefore, in 
many image processing problems, one can consider a mini- 
mization problem with any distribution if the statistical fea- 
tures of the underlying process are known, e.g. Poisson for a 
counting process, Rayleigh for ultrasound images, and so on. 

This paper is organized as follows: the region-based vari- 
ational framework for segmentation is introduced in Section 
121 Influence of the parameters estimator is investigated in Sec- 
tion [3] Our new theoretical results are detailed in Section |4l 
Experimental results are discussed in Section |5j We finally 
conclude and give some perspectives. 

2. PROBLEM STATEMENT 

Let U he a class of domains (open, regular bounded sets, i.e. 
C^) of M*^, and Hi an element of W of boundary dili. The im- 
age domain is noted flj. A region-based segmentation prob- 
lem aims at finding a partition of fi/ in n regions { i , . . , ri„ } 
of respective boundaries {dili, ..,dfln} that minimizes the 
following criterion: 

n 

s(f]i,..,r!„,r) = ^s„(r!,) + A£;6(r) (i) 

where Eri describes the homogeneity of the region 57^, and 
r = UiLi The energy term Ei, is a regularization term 
balanced with a positive real parameter A. 

2.1. Statistical energy 

Let us denote Y(x) a vector of random variables (VRV) at 
location x where x G 17/. A region Hi is considered to be ho- 
mogeneous according to the image feature y(x) e x C M'' if 



the associated VRV Y(x) follows a prespecified probability 
distribution with pdf p(Y(x), 6i) of parameters 6i in the re- 
gion Vli. The image feature may be for example the intensity 
of the region (a 3-dimensional vector for color regions d = 3, 
or a scalar for grey level regions d = 1). We then consider 
functions of the form: 



E„{ni)^ ( $(p(y(x),6»,)dx 

Jn, 



(2) 



with $ at least and Lebesgue integrable function. 

Since the above integrals are with respect to domains, and 
also that the parameters Oi are region-dependent, optimiza- 
tion (and hence derivation) of (|2]l is not straightforward. The 
derivation of such a functional is performed using domain 
derivation tools as in |4j. The following theorem from [4J 
will be useful for region-based terms derivation: 

Theorem 1 The Gateaux derivative of the functional 

J{fl) = J /(x, fl) dx in the direction ofV is the following: 

n 

<j'{n),v>^ J f,{^,n,v)d^-J /(x,i])(v.Ar)da(x) 

a an 

where N is the unit inward normal to dQ, da its area element 
and fs the shape derivative of f. 

2.2. Towards a geometrical PDE 

From the shape derivative, we can deduce the evolution equa- 
tion that will drive the active contour towards a minimum of 
the criterion. In order to fix ideas, let us consider the case of 
a partition of an image in two regions {Oi , 172}. The curve F 
stands for the interface between the two regions. The bound- 
ary energy term Ei, may be chosen as the curve length and 
derived classically using calculus of variation |7| or shape 
derivation tools |4|. 

Let us suppose that the shape derivative of each region fl^ 
may be written as follows: 

<EUm,V>^~ [ /,(x,l],)(V(x).Ar(x))da(x) 



(3) 

In this paper, we focus our attention on the computation 
of the shape derivative of the criterion (|2]). Let us note that 
when — — log(t) the function ^ is known as the log- 
likelihood score function, used to describe the homogeneity 
of a region according to a prescribed noise model. It has 
been used in |1, 2, 3|. Under some conditions, we are able 
to prove that minimizing this function with the ML estimator 
gives a simple evolution speed whose closed-form expression 
depends only on the pdf. However, one must be aware that 
using a more general function $ or different hyperparameter 
estimator will generally yield complicating additive terms in 
the shape derivative and hence in the evolution speed expres- 



3. INFLUENCE OF THE ESTIMATOR ON THE 
EVOLUTION EQUATION 

This section is devoted to illustrate two examples that support 
our previous claims on the influence of the parameters esti- 
mator More precisely, we state two derivative results in the 
case of the Rayleigh distribution with two different estimators 
(ML and moments). 

3.1. Shape derivative for the Rayleigh distribution 

We consider the Rayleigh distribution of scalar parameter 6*^: 



i t \ a \ y w ( -y(x)^ 



2ef 



(4) 



Classically, one can compute an estimate of the parameter 6i 
using the moment method. In this case the estimator is given 
by: 



where = Jq. dx and y(x)^ denotes the sample mean 
inside the region (li. 

Alternatively, one can also compute an estimate using the ML 
estimator given by: 



2\ni\ 



y(x)2dx 



(6) 



3.1.1. Shape derivative with the moment estimator 



Theorem 2 The Gateaux derivative, in the direction of V, 
of the functional Eri{Qi) = - J^^ log(p(y(x), 6iMo)dsi{yi) 
with p a Rayleigh distribution, is the following: 



<E'„{n,),V> = / (log (p(y(x),0, mo) 
+ ^(y(x),l],)(V- Ar)da(x) 



7i"y^(x), 



y(x). 



1 - 



y(x) 
y(x). 



where y'=(x), = J^^ y'=(x)dx. 



3.1.2. Shape derivative with the ML estimator 

Theorem 3 The Gateaux derivative, in the direction ofV, of 
the functional Eri{fli) = — Jq. '^og{p{y{x), 9iML)d'X- with p 
a Rayleigh distribution, is the following: 

<EU{n,),V>= f log (p(y(x),0, ml)) (V-N)da(x) 

One can then see that the parameter estimator has a clear 
impact on the evolution speed expression, as the additive term 
A(j{x) , Hi ) appears in the shape derivative when the moment 
estimator is used. 



4. GENERAL RESULTS FOR THE EXPONENTIAL 
FAMILY 

Let us now consider the shape derivative within the frame- 
work of muhi-parameter exponential family. This family in- 
cludes Poisson, Rayleigh, Gaussian... 

4.1. Exponential families 

4.1.1. Definition 

The multi -parameter exponential families are naturally indexed 
by a fc-real parameter vector and a fc-dimensional natural statis- 
tics T(Y). A simple example is the normal family when both 
the location and the scale parameters are unknown {k — 2). 

Definition 1 The family of distributions of a VRV Y {Vrj ■ 
T] £ £ C R''}, is said a A;-parameter canonical exponential 
family, ;/ there exist real-valued functions r/i , . . . , ry^ : O i— > 
R and A(r]) on £, and real-valued functions h,Ti, . . . ,Tk : 
R'* M, such the pdfp{y, ij) of the Vrj may be written: 

P{y,r]) = /i(y)exp[(r,,T(y)) - A(r7)],y e x C R'^ (7) 

where T = (Ti, . . . , Tk)^ is the natural sufficient statistic, 
T] = {rji, . . . , iik)^ and £ are the natural parameter vector 
and space. (77, T) denotes the scalar product. 

We draw the reader's attention to the fact than 77 is a function 
of 9 which is the parameter of interest in most applications. 
For example, for the Rayleigh distribution : rj — 

T(y) = y\ h{y) = ^ and A{rj) = - log(-2r,). 

4.1.2. Properties 

The following theorem establishes the conditions of strict con- 
vexity of A, and then those for VA to be 1-1 on £. This is 
very useful result for optimization (derivation) purposes: 

Theorem 4 Let V a full rank {i.e. Cow[r(Y)] is a positive- 
definite matrix) k-parameter canonical exponential family with 
natural sufficient statistic T(Y) and open natural parameter 
space £ [S]. 

(i) VA : f I— > 5 C R ii 1-1. The family may be uniquely 
parameterized by fJ,{r]) = £^(T(Y)) = VA(rj). 

( ii) The -log-likelihood function is a strictly convex function 
of r] on £. 

These results establish a 1-1 correspondence between r] 
wdE{T{Y)) such that: 

S3 11 = VAir]) = E{T{Y)) ^ £ 3 r] = i; {E{T{Y))) 

(8) 

holds uniquely with and continuous. 



4.2. General results 

In the sequel, for the sake of simplicity, we will invariably 
denote rj for the natural parameter and its finite sample esti- 
mate over the domain (without a slight abuse of notation, this 
should be r)). 

Theorem 5 The Gateaux derivative, in the direction ofV, of 
the functional Eri{fli) = /j^, $(p(y(x), rj(f2i)))c?a(x) where 
p{.) belongs to the multi-parameter exponential family and is 
expressed by and rj the natural hyperparameter vector, is: 

< E'„{n,),V >=- f <P{p{y)){V • N)da(x) 
Jdn 

+ f p(y)<i>'(p(y))(Vvr/,T(y)-VA(j7))dx (9) 

Jn, 

with Vv^y the Gateaux derivative oft] in the direction ofW, 
and (x, y) the scalar product of vectors x and y. 

In a finite sample setting, when using the ML estimator, 
we can replace VA(r7) by T(Y) (the 1st order sample mo- 
ment of T(Y)). Thus, when using the -log-likelihood func- 
tion, the second term becomes equal to J^^ (Vv^7, T(y) — 

T(Y)))(ix, and hence vanishes. The following corollary fol- 
lows: 

Corollary 1 The Gateaux derivative, in the direction ofY, of 
the functional Eri{VLi) = - /^^ . log(p(y(x), r5ML(^^j))'^aW 
when 77ML is the ML estimate, is the following: 

<E'„m,Y>= f (log(p(y(x),77ML(r!,)))(V.Ar)da(x) 

Jani 

This provides an alternative proof to the result of IT] [51 . 
Nonetheless, we here point out that, in the work of IT] [3], the 
role of the parameters estimator was not elucidated. 

5. EXPERIMENTAL RESULTS 

This section presents some experimental results on noisy im- 
ages. The initial noise-free image is shown in FigHI For four 
different Battacharya distances (BD), we have systematically 
corrupted this image with two types of noise: Poisson and 
Rayleigh. The Battacharya distance is used as a measure of 
"contrast" between objects and background. It is defined as : 

T^{pf{y):Po{y)) = - log / Jpf{y)po{y)dy 

For each combination of BD value and noise type, 50 noisy 
images were generated. Each noisy image was then segmented 
using four different energy functionals, namely Chan-Vese 
||9l , and our method with -log-likelihood and ML estimator 
with three assumed noise models: Gaussian, Rayleigh and 
Poisson. For each segmented image with each method at each 




Fig. 1. Initial noise-free image. 




Fig. 2. FPF and TPF as a function of BD. 



BD value, the average false positive fraction (FPF) and true 
positive fraction (TPF), over the 50 simulations were com- 
puted. The bottomline of these experiments is to show that 
using the appropriate noise model will yield the best per- 
formance in terms of compromise between specificity (over- 
segmentation as revealed by the FPF) and sensitivity (under- 
segmentation as revealed by the TPF). 

Fig m depicts the average FPF (left) and TPF (right) as a 
function of the BD for Poisson ((a)-(b)) and Rayleigh ((c)- 
(d)) noises. As expected, the FPF exhibits a decreasing ten- 
dency as the BD increases, while the TPF increases with BD, 
which is intuitively acceptable. More interestingly, the best 
performance in terms of compromise between FPF and TPF 
is reached when the contaminating noise and the noise model 
in the functional are the same. This behaviour is more salient 
at low BD values, i.e. high noise level. One can also point 
out that the Chan-Vese functional is very conservative at the 
price of less sensitivity. Clearly this method under-segments 
the objects. 



6. CONCLUSION 

In this work, we proposed a novel statistical region-based ac- 
tive contours method, where the region descriptor is written 
as a function $ of some pdf belonging to the exponential fam- 
ily. The case of the likelihood score, which describes the ho- 
mogeneity of the regions, is obtained as a special case of our 
setting. We shedded light on the influence of the noise param- 
eters estimator, and on the adequacy between the noise in the 
observations and the one incorporated in the active contours 
functional. Our ongoing work is directed towards extending 
our approach to n regions by adapting the multiphase method 
|9|. Furthermore, encouraged by our preliminary application 
results, a deeper experimental work is currently carried out to 
validate our approach on real data sets. 
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